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Abstract 

The phase transition for a finite volume system that incorporates the Polyakov loops and main- 
tains the colorless state is explored using the Polyakov-loop extended Nambu-Jona-Lasinio (PNJL) 
model. The order parameter for Polyakov loops is demonstrated to signal the appearance of a tran- 
^vq ' sition for SU{3)c analogous to Gross-Witten (GW-) phase transition instead of the deconfinement 

phase transition to quark-gluon plasma. The asymptotic restoration of Polyakov loops is conjec- 
tured to be a threshold production for meta-stable Hagedorn (or semi-QGP) states and this does 
^ ' not imply a direct deconfinement phase transition. In this context, the GW-like point is the point 

a: 

where the colorless states switches from the low-lying hadronic states to the meta-stable high-lying 
►^ , Hagedorn states. The chiral phase transition takes place within an extended GW-like point de- 

Qs^ I pending on the fireball's size. The deconfinement phase transition is determined by Hagedorn's 

^ ■ 

1/^ ■ temperature above GW-like temperature. 
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I. INTRODUCTION 

Recently, Fukushima [l| has extended the Nambu-Jona-Lasinio (NJL) model to include 
Polyakov loops, namely, $ and $ and the a-chiral field. Fukushima's approach is known as 
Polyakov extended Nambu-Jona-Lasinio (PNJL) model. The Polyakov loops are related to 
the imposition of the Gauss' law where the trivial vacuum is a minimum of the free energy 
and the formation of stable colorless QG-droplet [2]. The PNJL model has been widely 
adopted to study the phase transition diagram. Furthermore, it has been extended to 
investigate the phase transition diagram with various phenomenological effective Polyakov 
and gluon potentials as well as various extensions to include other NJL's fields such as 
the isospin scalar and vector fields as well as the color superconductivity [jj, |3|4l3|. The 
hybrid description of PNJL with two and three fiavors has been studied extensively with 
various modifications of the effective Polyakov and gluon potential and the results have been 
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compared with lattice QCD data |^. The Polyakov potential as a function of Polyakov loops 
(or equivalently the VanderMonde potential as a function of fundamental eigenvalues of the 
Gauss-law) is originated from the invariance Haar measure of SU{Nc) in order to project 
;he colorless state of the quark and gluon fireball. The comparison with lattice calculations 



i, y, y, Is 



hints that the effective Polyakov potential could by modified by temperature. 
Furthermore, it has been suggested that modifying the invariance Haar measure's exponent 
somehow mutates the Hagedorn's internal structure [ij]. Furthermore, the bag's volume 
fiuctuation can modify the effective Polyakov potential. Nevertheless, the modification of the 
effective Polyakov potential in the medium and its impacts in the phase transition diagram 
will be considered in a future work. The internal structure of the quark-gluon (QG) has 
been suggested to be crucial to the tri-critical point (see for instance Ref. 15| and reference 
therein). This has a significant impact in the recent research to explore the width of phase 
transition and intermediate processes such as Hagedorn states, quarkyonic matter and semi- 
classical QGP phases and color-fiavor superconductor matter. There exist various reviews 
discussing the QG-blob's internal color structure. Brezin, Itzykson, Parisi and Zuber studied 
;he planar approximation to field theory through the limit of a large internal symmetry group 



16|. This procedure is known as the matrix saddle point method. Gross and Witten [17 1 
using the matrix method have discovered a possible transition from a specific phase with 
strong coupling to another phase with weak coupling in the large Nc limit (i.e. Nc ^ oo 



but a finite g'^ Nc) of Wilson lattice gauge theory. For technical reasons, the spectral density 
method which has been developed by Brezin et. al. |l6|, ll7| depends basically on the 
large N^ limit and it is not permissible for technical reasons to extend the same analyses 
using the spectral density for finite number of colors. The GW-like point sticks in one's 
mind for N^. ^ oo and remains obsolete for N^ = ?> (i.e. the QCD). The GW-like phase 
transition for finite N^ is not expected to have the same characteristic behavior to that one 
in the limit N^ — )■ C)0. In order to search for a mechanism analogous to GW-transition in 
QCD, it is important to extend the analysis using the (non-Gaussian-) stationary points 
method in the strong coupling limit in the context of the Polyakov loop parameterization 
as done by in Ref jl| and the references therein on one hand and the (Gaussian-) saddle 
points approximation in the weak coupling limit as done by Elze, Greiner and Rafelski 
and others 
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on the other hand. The interpolation between the asymptotic non- 
Gaussian stationary points approximation's solution for the strong coupling limit and the 
asymptotic Gaussian saddle points approximation's solution for the weak coupling limit 
is not fully understood in QCD and the corresponding mechanism is analogous to GW- 
transition. Furthermore, GW-like transition sounds to take place over the interpolation 
range between two asymptotic solutions (i.e. over an extended interval) rather than a single 
defiection point. Furthermore, Elze, Greiner and Rafelski have pointed out that the non- 
perturbative effect of the colorless state leads to a gradual freezing of internal degrees of 



freedom 20|, |21|. This mechanism could explain the emergence of QG liquid droplet(s) or 
equivalent forms such as Hagedorn states, quarkyonic droplets etc. It should be stressed 
that GW-like transition is not a confinement /deconfinement phase transition, but instead 
is the production threshold of (meta-) Hagedorn states in hadronic matter. The Hagedorn 
states emerge as gas of bags. Therefore, there is a possibility for a new form of matter that 
can be formed in a narrow range above GW-like point and below Hagedorn's temperature. 
This form of matter emerges as a gas/liquid of bags and these bags expand and grow up 
gradually. When Hagedorn's temperature is reached, the system undergoes a deconfinement 
phase transition to QGP. 

The outline of the present paper is as follows: In Sec. 11, we review Polyakov loops 
without chiral field and demonstrate the interpolation between the low-lying and high-lying 
energy solutions and a possible transition that is analogous to GW-transition. In Sec. HI, 
the treatment is extended to include the cr-chiral field in the context of PNJL model and 



demonstrate the emergence of an extended GW-like point. The connection between GW-hke 
point and production of Hagedorn states is discussed in Sect. IV. Finally, the conclusion is 
presented in Sec. V. 

II. A SIMPLE CANONICAL ENSEMBLE WITH POLYAKOV LOOPS 

The grand potential for the quark and the anti-quark is given by 
V 
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where e^ {p) = ^/p^ + m^, (2 J + 1) = 2 is the spin degeneracy, V is the quark and gluon 
blob's volume and fig is the flavor chemical potential while di are the imaginary color chemical 
potentials or fundamental gauge fields of the SU{Nc) group's fundamental representation 
(i.e. Gauss-law's eigenvalues on the thermal excitations). When no chiral fields are involved 
in the calculation, mg is reduced to the current mass (for only the sake of simplicity, it can 
be assumed massless for light flavors). The first term in the square bracket that appears on 
the right hand side of Eq.([l]) is temperature independent. It diverges at zero temperature 
and is a non-re-normalizable term. It can be regulated in the standard way by introducing 
UV-cutofF for the momentum integration. In the standard a-model, that term is trivially 
dropped as far it can be absorbed by the nonlinear cr-potential but the this is not the case 
in Nambu-Jona-Lasinio model (NJL) where the first term is regularized and retained in the 
calculation. After a simple algebraic manipulation, Eq.(Il]) becomes 
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where A is UV-cutoff that regularizing the divergent term over the niomentum integration. 
The UV-cutoff for momentum integration is taken A = 631.5 MeV in the present calculations. 
The Polyakov-loop triality parameters $ and $ are defined, respectively, as follows 

$ = -^ [e-''' + e-''' + e-''-'] , (3) 

where 0^ = —9i — 02 for SU{3)c and the fundamental gauge fields 9i, i = 1,2,3 are subjected 
to the periodicity condition over the interval — vr < 6i < tt. Eq.([2]) can be written as follows 
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In the case of massless flavors and /x^ = and in the terms of fundamental gauge fields, 
Eq.dl]) reads 
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The canonical ensemble for the quark and anti-quark becomes 

where fi^g (/3, y; $, $) is given by Eq.dl]). Fortunately, the qq grand potential becomes a 
real one when /i^ = 0. The partition function for the gluons can be calculated in a similar 
manner. The phenomenological gluon potential parameterized in the terms of Polyakov 
loops has been adopted recently in the literature |l|. The general choice is given by 

lfi,(/3,K;$,$) = -iilog,Z,(/3,l-;$,$), 



where Polyakov term $$ can be written in the terms of fundamental gauge fields as follows 
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where Oq, fli and 02 are phenomenological parameters. The phenomenological gluon partition 
function is usually adopted in the term of Polyakov loop approach in order to replace the 
standard gluon partition function that is given by 
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where e^ 



|p|. Eq. (lTO|) is evaluated explicitly as follows 
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In the standard treatment the gluons are treated as the adjoint interaction particles of the 
SU{Nc) symmetry group. It should be noted here that the SU^NcYs adjoint eigenvalues (i.e. 
adjoint gauge fields), namely, 0" are calculated from the nested commutation relations for the 
fundamental eigenvalues, namely, 6i of the Lie algebra. The adjoint eigenvalues are related to 
fundamental eigenvalues by the relation 0° = {6i — Oj). This relation diagonalizes the adjoint 
representation and subsequently commutes with the Hamiltonian. In order to understand 
the origin of the gluon's phenomenological potential, Eq.l ITOl) can be approximated and 
simplified in order to be evaluated using Polyakov loop variables in the following systematic 
way 
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where the factor 2 that appears on the right hand side conies from the spin degener 



acy. The coefficients C„ are functions of Polyakov loops 32|. When the Polyakov loops 



vanish ($, $ — )■ 0), the gluon grand potential is reduced to \im.^-^_^^yVlg (/3, V; $, $) = 
^ ^ ^ " ^ while in the case of Polyakov loop restoration ($, $ — > 1), it is reduced 
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to Mm^-^^^^ilg (/3,\/;$,$) = —ji "45 — • This implies that Q,g is reduced by factor 
l/N^ when $, $ are changed from 1 to 0. It can be parameterized to yflg (/3, V;^,Wj ^ 

~W ^ N^ " 4g — . In order to simplify the calculation drastically, the gluon grand 

potential is simplified to a phenomenological potential such as that one given in Eq.([7]) as 
follows 
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In the calculation of the phase transition from the low-lying energy excitations to the high- 
lying ones but below the deconfinement phase transition, it is adequate to use the potential 
that is given by Eq.f lT^ . 

The canonical ensemble for a finite volume quark and gluon blob in the Hilbert space is 
given by the Fock product of quark and antiquark partition function and the gluon partition 
function as follows 

Zqqg (/3, V] 9i, 62) = Zqqg {(3, V\ 9l, 62, 03 = —&! " ^2) , 

= Z,5,(/3,V^; $(^1,^1), ^(^1,^1)), 

= Z,g(/3,V^;$,$)xZ,(/3,V;$,$). (15) 

This implies that the grand canonical ensemble is reduced to 

Z,53(/3,V;<I>,$) = exp(-/3fi,g,(/3,V; $,$)), 

= exp(-/3fi,5,(/3,V;^i,^2)), (16) 

where 

f),5,(/3,V^;$,$)=3?ef],5(/3,V;<l>,$)+fi,(/3,V;<l>,$). (17) 
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The colorless state for the quark and gluon blob is ensured by projecting the color singlet 
state in the following way 
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The VanderMonde potential is stemmed from the invariance Haar measure of the group 
integration and is defined by 
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where the parameter Ggym depends basically on the group's symmetry. It is reduced to 
Gsym = 2 for SU{Nc). In the lattice modehng, the number of states for the VanderMonde 
potential |2| |33[ is introduced by / d'^xd''' (0) as follows 
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where a is the lattice size. Hence, the VanderMonde potential is regulated [2,] as follows 
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2 sm ' 



(22) 



A finite bag with volume at the same size order of the lattice l^ ~ a^ ~ fm and ^^ ~ 1 
is considered in the present work. The regulation 7^63 = ^ will be considered elsewhere. 
However, the term ^\^ in MIT bag model is related to the volume fluctuation for a bag with 
an extended surface. Nonetheless, the VanderMonde potential regulation is essential for a 
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system with infinite volume |2|. Hereinafter, tlie number of states is considered -\V = 1 for 
VanderMonde potential in colorless quark and gluon bag. 

In order to consider Polyakov loops parameterization, it is useful to perform the variable 
transformation from fundamental gauge fields, namely, {61,62,63) with 6^ = —61 — 62 to 
Polyakov loop variables, namely, ($,$). In the context of SU{3)c, the invariance Haar 
measure is furnished by 

Ms) = TTTTT^ / d6, f d62 J] 2sin "" 



A^! (27r)2 
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The invariance Haar measure can be transformed and written in the terms of Polyakov loop 
variables $ and $. The transformation of the square root of the invariance Haar measure 
from the variable set {6i} to $ and $ leads to 

'6,-6- 
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The integration over 61 and 62 is transformed to Polyakov loop variables $ and $ as follows 
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Hence, the invariance Haar measure that is given by Eq. (ll9p becomes 
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N !(27r )^c-i ^'-'^ SU{?>)c- The subscript notation, namely, C that appears 
under the integral indicates the integration is over a complex plane domain. The complex 
domain for $ and $ is the three pointed star with a radius 1. The complex domain for 
Polyakov loops complicates the situation when the non-Gaussian stationary point method 
fails and the Gaussian saddle point procedure turns to be essential. The invariance Haar 
measure with Polyakov loops parameterization can be represented as an effective Polyakov 
VanderMonde (PVdM) potential. The colorless canonical ensemble with an effective PVdM 



potential in SU{3)c group representation reads 

Zcolorless (/?, V) = N^^^^ J d$ J rf$ [l - 6$$ + 4 (^^ + $') - 3 ($$)' 

X Z,5,(/3,y ;$,$), 
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where PVdM potential in SU{3)c representation is given by 



VpydM (/?;$,$) = -^^log, (l-6$$ + 4(<l>3 + $')-3(<l><l>)'), (28) 



and 



n,^, (/3, V;^,W} = -i log, [Z,^, (/3, V; $, $)] . (29) 

The phenomenological PVdM potential can be introduced by adding a phenomenological 
pre-factor parameter, namely, aph, in front of the logarithm as follows 

Vpv,m(/3;$,$) = -iiap/.log, (l-6<l>$ + 4(<l>=^ + $')-3(<l>$)'). (30) 

This phenomenological parameter, namely, aph modifies the underlying internal symmetry 
of Hagedorn states and in some scenarios this could break the internal symmetry of the 
QG-bags but not the global symmetry of the system [ij]. The phenomenological PVdM 
potential and the variation of the phenomenological parameter, aph, will be considered in 
another work. The integral that is given by Eg. (1271) is evaluated using the non-Gaussian 
stationary points method over the complex plane. The Polyakov's stationary points, namely, 
$ = $0 and $ = $0 are evaluated by extremizing the exponent term. The stationary points 
$0 and $0 are calculated as follows 



^ ^ [VpvaM (/3; $, $) + fi,5. (/3> V; $, $)] 1^^. ^^^^ = 0, 



' ' [Vpv.M(/3;$,$)+f^,,,(/3,V^;$,«l')]L^^„^^^ =0. 



The $'s extremum is determined by the following constraint. 
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The same thing can be done for $. As far as the nuclear matter environment remains in the 
circumstance that Polyakov's stationary points are located in the region $o < 1 and $o < 1 
(i.e. non-Gaussian stationary points) and below the threshold of GW-like phase transition, 
then the canonical ensemble which is given by Eq. (p7|) . is evaluated as follows 

Zilriess (/3, V) = exp (-/3 [Vpv,M (/3; $o, $o) + ^,-,, (/3, V; $o, $o)]) • (33) 

The pre-factor constant, namely A^Haar' ^^^^ appears in Eq. (!27|) is dropped in order to 
normalize the partition function. The solution that is given by Eq.f l33|) is assigned as the 
low-lying energy solution (I). This solution is the asymptotic solution below the threshold of 
GW-like phase transition point. The validity of the low- lying energy solution (I) is satisfied 
as far Polyakov's non-Gaussian stationary points remain in the energy domain $o << 1 and 
$0 << 1 (i-e. far away from Polyakov triality restoration point). In the case that fig = 0, 
then the equations' set given by Eq. fl32l) becomes symmetry over Polyakov loop variables and 
this leads to equal and real stationary points for $ and $. However, whenever |$|o — ?■ 1~, 
then the effective PVdM potential, namely VpvdM (/?; ^, *&), develops a virtual logarithmic 
divergence. Therefore, the logarithmic divergence of the effective PVdM potential spoils 
badly Polyakov's non-Gaussian stationary points procedure and leads to virtual singularity 
for the effective grand potential of the system. Evidently, this virtual singularity deforms 
the low-lying energy solution Z^J^^^^^^ (/3, V). This kind of behavior indicates modification 
in the analytic behavior of the canonical ensemble and another analytical solution, namely 
the solution (II) emerges in the system. The change in the analytical solution is the beneath 
mechanism of GW-like phase transition even for finite number of colors (i.e. Nc = 3). 
The GW-like point for A^^ = 3 may play a significant role in the deconfinement phase 
transition diagram in nuclear physics as far it is not a confinement /deconfinement point. At 
the onset of GW-like phase transition, the non-Gaussian stationary points turn to behave 
as Gaussian saddle points that oscillate harmonically around the stationary points. This 
mechanism reflects the modification in the analytical behavior from the asymptotic solution 
(I) to the solution (II) when the temperature reaches GW-like point. Therefore, when GW- 
like threshold is reached, the canonical ensemble modifies its characteristic behavior from the 
low-lying energy solution ^^oL/e^s (Z^' ^) ^o the high-lying energy solution Z^JJli^^^{P,V). 
The second solution implies a possible production of Hagedorn states. When the asymptotic 
high-lying energy solution ( i.e. solution (II) ) is reached, it becomes more suitable to 
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write the invariance Haar measure in the terms of fundamental gauge fields (i.e. 6i, 62, 
^3 = —9i — 62) rather than Polyakov loop variables (i.e. $, $). Therefore, at the threshold 
of GW-like phase transition, the colorless canonical ensemble ( i.e. solution (II) ) is reduced 
to 
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Furthermore, in the case of massless flavors and zero flavor chemical potential (i.e. /Xg = 0), 
the canonical ensemble (II) for the high lying energy solution is simplifled to 



(37) 



In the present model, the order parameter(s) of GW-like phase transition is (are) tem- 
perature (and/or flavor chemical potentials). At that critical point T = T^., the low- lying 
and high-lying energy solutions match each other. In this case, the low-lying energy solution 
is extrapolated to the high-lying energy solution at the threshold of GW-like point. The 
critical value of T^ = l//3c is determined by the continuity condition 
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Below GW-like point, the non-Gaussian stationary point of the asymptotic solution (I) is 
limited to |$| < 1. At the threshold of GW-like phase transition, the non-Gaussian station- 
ary points turn to be Gaussian saddle points that oscillate in the neighborhood of the center 
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of the symmetry group. Hence, the solution (II) turns to the asymptotic solution above GW- 
like point. The exact solution of Eq.f llSp is obtained by evaluating the integration over the 
invariance Haar measure numerically. It is found that the asymptotic solution (I) matches 
the exact numerical solution below GW-like point while the solution (II) matches the exact 
one above GW-like point. Solutions (I) and (II) intersect each other in the neighborhood of 
GW-like point. The Gaussian saddle points procedure is better understood in the terms of 
fundamental gauge field variables 9i rather than Polyakov loop variables $ and $. Further- 
more, the deconfinement phase transition takes place when the high-lying energy states of 
quark-gluon bags become unstable and in this case the Hagedorn matter undergoes phase 
transition to quark-gluon plasma at Hagedorn's temperature. It is worth to note here that 
below Hagedorn's temperature, the high-lying energy quark-gluon bag acts as quark-gluon 
fluid (or semi-QGP) droplets as far the constituent quarks and gluons remain within the 
range of the effective VanderMonde potential interaction. 

III. THE EXTENSION TO POLYAKOV-NAMBU-JONA-LASINIO MODEL 

The conventional NJL Lagrangian density reads 

Cnjl = q [t^d^ - m,] g + ^G [{qqf + (qt^.Tq)'] , (39) 

where ruq is the quark's current mass and G is the NJL coupling constant. The constant 
G is adjusted in order to fit the nuclear phenomenology. The current mass and coupling 
constant for light flavors are taken mg = 5 MeV and G = 10.992, respectively. The quark 
and antiquark grand potential density in the presence of the effective chiral field is furnished 

by 

where Eq^{p) = ^Jj? + Mq{aY and Mq{a) = m^ — Ga. The parameters m^, a, $, G and 
liq are quark's current mass, scalar field, Polyakov loop parameter, scalar coupling constant 
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and constituent quark's chemical potential, respectively. The a scalar mean field indicates 
the condensate a =< qq>. The first term on the right hand side of Eg. (1401) is temperature 
independent and is regularized as follows, 
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where the UV-cutoff A = 631.5 MeV. The scalar a-chiral field in the context of PNJL model 
is considered self-consistently. In the terms of fundamental gauge fields {6i, 02, 9^ = —61 — 62) 
rather than Polyakov loops ($,$), Eg. (140 p is reduced to 

fi,5(/3,y;a,<l>,$) nq-q{(3,V-a,'^{6^,62)M6^,62)) 
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(],g(/3,y;a,gi,g2 
V 



(42) 



where 



V 



N 



f rA 



2iVcE 



d\p\\p\\ 
27r2 '"'^'^ 



q k=l 
+ 



1 + cos(^fc) e/^[^«-(P^-^9] 



(27r2) 3Eg^ (p) Vl + 2 cos(^fc) e/^[-^9-(p)-/^«l + eS/^I-^^-CpO-/^^] 
1 + cos(^fc) e''[-^^-(P^+^9 



1 + 2 cos(6'yt) e/^['^9<T{R)+M,] + (,2p[Eg^{p)+fi. 

q k=l 



smWk 



(27r2) 3Eg„ (p) V 1 + 2 cos(^fc) e/3[-E^.-(pO-M,] + emE,.{p)-f^,] 
sin (6k) 



g/3[£,,(p)+A.g] 



1 + 2 COS {6k) e^[^9-(P^+^9] + e2/3[S9<T(p)+M9] 
Therefore, in the case fig = 0, Eg. (j43l) becomes real and is simplified to 



(43) 



ng-gif3,V;a,6,,62) 
V 

q k=l 



N 



f /.A 



27V. 5:/ 

„ "'0 



d\p\ bT 



d\p\\p\\ (-^ 



l + cos(^fc)e^^«-(P^ 



1 + 2 cos(^fc) e^^'^'^Pi + e2/5-^.-(p) / ' 



(44) 



(27r2) 3^,, (p 

Subseguently, the canonical ensemble for guarks and gluons in the context of the PNJL 
model becomes 



Z,5, (/3, V- a, $, $) = exp [-/^f],^, (/3, V- a, $, $)] . 



(45) 
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The total grand potential for cliiral quarks and gluons reads 

3f?efi,5,(/3,y;cr,$,$) = fi,^ (/3, a; $, $) + fi, (/3, V; $, $) + Vf/ (a) , (46) 

where V is the system's volume. The chiral quark and antiquark grand potential 
f2qg(/3, V; a, $, $) is determined by Eq. fllUl) while the gluon grand potential ^g{(3, V; $, $) is 
determined from Eq.f ll3l) . The effective chiral potential is given by 

U{a)=^-Ga'. (47) 

The canonical ensemble for the colorless quark and gluon blob is determined by projecting 
the color-singlet state in the following way, 

Zcoioriess (/?, V ] a) = / dfi (g) Zggg (/3, V; a, $, $) , 

= jdixig)Z,-,gi(3,V;a,9,,e2). (48) 

It is possible to write Eq. (HHj) in the terms of Polyakov loop variables $ and $ as follows 

Zcoioriess (/?, V; a) = N^^^^ [ d^ [ d^ZpNjL (/?, V' a, $, $) , (49) 

Jc Jc 

where the subscript C indicates the integration over the three pointed star boundary in the 

complex plane and 

ZPNJL (/3, V; a, $, $) = exp {-/3 [VpvdM (/?; $, $) + ^,5. (/?, V; a, $,$)]) . (50) 

The double integrations over Polyakov loop variables, namely, $, $ is evaluated using the 
non-Gaussian stationary points method for the low-lying energy limit. Subsequently using 
the non-Gaussian stationary points method, Eq. (l49|) is reduced to the low- lying energy 
solution as follows, 

Zilriess (/3, V; a) = Ziil,^,, (/3, V; a, $0, $0) , 
= ZpNjL{l3,V;a,%,%) , 
= exp (-/3 [VpvdM (/3; $0, $0) + ^g-gg (/3, V; a, $0,^0)]) • (51) 

The pre-factor A^fja^j- is eliminated in order to guarantee the normalization of the non- 
Gaussian stationary point method and it does not affect the calculation. The Polyakov 
loops' non-Gaussian stationary points, namely, $0 and $0 are determined by extremizing 
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the exponent which appears on the right hand side of Eg. (1501) with respect to $ and $ in 
the following way, 



^ [VpvdM (/3; $, $) + n,-,g (/3, V- a, $, $)] 



a* 



iJ>=$Q,cJ>=<J)o ' 

0. 



(52) 



l'I'='I>0,*='J'0 

Furthermore, the a-chiral field stationary point, namely Uo, is determined by extremizing 
the exponent which appears on the right hand side of Eq. ( ISTl) with respect to the scalar field 
a. Since Polyakov VanderMonde potential, namely, VpvdM (/?; "^j '^) does not depend on a, 
the variation of ^qqg (/?, V; a, $0) "^o) with respect to a mean field leads to 



a 



1 d / l]gg(/3,l^;a,$o,$o) 



(53) 



where 



d / l]gg(/3,V-;a,$,'l>) 
'9a I l^ 



Aff 



E ^''■ 



■,('')^a/,(ct: 



CT=(TO 



2A/', 






6 



-6 



rf|pi IpI^ 



($ + 2^e-'^^'"'^P^) e-^'^""^^ + e"^'^^" 



■(p) 



27r2 £;g^(p) [1 + 3 ($ + <|)e-/3^9-{p)) e-^^i-^P^ + e-3/3^''-(p)] 



rf|p1 IpI^ 



($ + 2$e~^-^«'''^^) Q~l3Ega{p) _|_ ^-3l3Eg„{p) 



27r2 E^^(p) [1 + 3 ($ + $6-/3^5,. (p)) e-/3^«-(p) + e-3/3E,-(p)] 



(.54) 



Moreover, when Eg. (15^ is written in the terms of fundamental gauge fields {61,02), it is 
reduced to 

Nf 



d_[n^-q{l3,V- 0,6^,62, 
da 



V 



E (M,(a)^M,(a)) 



2Nr, 



d\p\ \p\' 



Nc 



-^E 



d\p\ \p\' 



2r^^E,M 
l + e^^^-(P^cos(^fc) 



^ io 27r2 Eg,(p) 1 + 2e/^-^«-(P^ cos (4) + e2/3-^«-(pT 



(55) 



The first term inside the square bracket on the right hand side of Eq. (15^ and/or Eq. (155|) is 
temperature independent. Its explicit expression reads 



2n^Eq,{p) 47r2 
1 



KjMJaY + A.^ 



An 



^M^yio^ 



A 



M,(a) 



A 



M,(0r 




(56) 
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Therefore, the partition function (i.e. the canonical ensemble) for the colorless quark and 
gluon bag with a-chiral field reads 

^colorless (Z^' ^) = ^colorless (/^' ^5 ^o) , 

= ^coLe.s(AV^;^0,$0,$0), 

= ZpNjL{P,V;ao,<l>o,%), (57) 

where the values of $0? ^o and o"o are the stationary points and they are calculated by 
Eqs.( !52|) and ( !53l) . respectively. Evidently, when the temperature approaches the critical one 
(i.e. GW-like point), the non-Gaussian stationary point method fails due to the logarithmic 
divergence of Polyakov VanderMonde potential, namely, VpvdM (/?; "^7 ^\ The logarithmic 
divergence of Polyakov VanderMonde potential indicates a collapse of the non-Gaussian 
stationary point method for the low-lying energy solution and the emergence of the high- 
lying energy solution where Polyakov loops' stationary points switch to become Gaussian 
saddle points that oscillate harmonically around the non-Gaussian stationary points. The 
Gaussian saddle point procedure is understood in the context of fundamental gauge fields 
9i and 62 much better than in the frame work of Polyakov loops $ and $. 

On the other hand, the partition function for the asymptotic high-lying energy solution 
(i.e. solution II) for colorless quark and gluon blob in the context of PNJL model reads 



Z^rless (/3, ^5 a) = / ^/^ (g) Zf^a {P^ V-^ ^, ^^ ^) 



1 1 



dOi / d^2 n (^* - ^i)' e-^^'^''('^'^'"'^i'^^), (58) 
-00 J —00 • ,• 



where ^3 = —Q\ — ^2- The quark and gluon grand potential which appears in the exponent 
in Eq.f l58|) is given by 

fi,5,(/3,\/;a,0i,02)) = l^,^3(/3,l^;a,<l>(^i,e2),$(^i,^2)), 

= fi,g (/3, V- a, ^1, ^2) + % (/3, V\ a, ^1, ^2) + y U{a), (59) 

where V is the bag's volume and Polyakov loop parameters $ and $ are written explicitly 
as functions of the fundamental gauge fields 6'i, 62 and 9^. The effective scalar potential, 
namely U{a), is given by Eq. (l47p . Since the non-Gaussian stationary points for the low- 
lying energy solution (I) are reduced to Gaussian saddle points for the asymptotic high-lying 
energy solution (II), it becomes essential to compute the quadratic expansion of the grand 
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potential around the Gaussian saddle points in order to evaluate the integral that is given in 
Eq.f l58|) more appropriately. This can be done much easier in the framework of fundamental 
gauge fields rather than Polyakov loops. The Gaussian saddle points of the fundamental 
gauge fields accumulate at the origin and fortunately this behavior simplifies the calculation 
drastically. The quadratic Taylor expansion of the grand potential density for quarks and 
anti-quarks is reduced to 

-n,-, (/3, V- a, ^1, e,) = -ng) (/3, V; a) + --fiJJ (/3, V- a) J^ ^^ (60) 

i 
th 

where the term reads 



:n^S{P,V-a) = -2N^J2l (^^-(^"^ 



N 



/■ d\p\ \p\^ 



_4NT r^ -^ - f61) 

^ ^ J 27r2 3Eg^{p) [e^E.^iP^ + 1] ' ^°^ 



while the quadratic term is given by 

Again, the gluonic grand potential for the low-lying energy colorless quark-gluon bags is 
assumed to be adjusted by the phenomenology as done in SecjIT] (see for instance Eq.( IT3|) ). 
This class of the phenomenological gluon potential is inspired from lattice calculations and 
has been recently adopted widely in the literature (for instance see [j]). The quadratic 
Taylor expansion of the gluon grand potential which is given by Eq.(TT3l) is approximated to 

^fi. (/3, V) = IfiW (/3, V) + ^ifif (/3, V) Y, {0. - O^f ■ (63) 

th 
The term reads 

i«f(ftV)^-i(|). (64) 

while the quadratic term is reduced to 
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where Ug = — ^^ — -. After evaluating the Gaussian integration over the fundamental 
gauge fields, the partition function for the high-lying energy solution (i.e. solution II) is 



approximated to 



^colorless (/^' ^) " ^colorless (/^' ^5 '^o) 



(66) 



where 



^(11) 



colorless 



W,V;a) 



)(0) 



(0) 



UZ,' n\ exp -/3 n%> (/3, V; a) + Q^ {f3, V) + V U{a) 



/iVe(27r)i(^=-i) 



r-H67) 



/3 



)(2) 



(2) 



n]i>{P,V;a) + 2N,n);>{P,V) 



Furthermore, cr-mean field (i.e. ao), above the threshold of GW-like phase transition, is 
determined by calculating a-stationary point in the following way. 



d n%> (/3, V; a) 



da 



V 



0, 



(68) 



a=cro 



where 



C^*(/3,V^;^) = fi;y(/3,V^;a) + fif (/3,r) + Vf/(a). 



(69) 



Under the assumption of the stationary point method, the extremization procedure is per- 
formed for the exponent term that appears in Eq. (1671) . The extremization of Eq. (1681) leads 
to 



Ga 



d_ ( n'S (/3, V; a) 
da 



V 



(70) 



where 



d_ ( n^' (/3, V; a) 
da 



V 



Nf 



2N^Y1 



dMJa 



d\p\ p^ 



27r2Eg,(p)(e/3^«-(p^ + l)_ 



x(M>. ^^ 



(71) 



In order to calculate other thermodynamics quantities, the derivative of the partition func- 
tion with respect to X is reduced to 






loge 



Ziilliess iP, V; a) 



d 



d 



d 



g^^ (/3, V; a) + ^Qf (/3, V) + ^ [V U{a)] 



+ '"-0^{^^°^^(^ 



n^^iP,V;a) + 2N,n'~'\p,V) 



(72) 
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where X is a thermodynamic ensemble such as /3 and V. For instance, from Eq. (172|) . the 
grand potential for colorless quark and gluon blob reads 



^C^(/5^^;-) 






loge 



^i?(/3,^;^) 



V 



ziilLss i(3, V- a) 



^mm^ui.) 



a 



V 



V 



(73) 



Hence, if the a-chiral field is not restored below GW-like point, then it will be a discon- 
tinuity (i.e. at least of a higher order discontinuity) from ctq = o"o to ctq = ctq i^ ^^^ 
neighborhood of GW-like point because the value of o"o below GW-like point is determined 
by Eq.(l53l) while ctq above GW-like point is determined by Eq.(l68l) or Eq.(I70l). However, 
the extensive numerical calculations show that the chiral symmetry restoration usually oc- 
curs in the neighborhood of GW-like point and before the ultimate point of the extended 
GW-range is reached. This indicates that the high-lying energy solution is chirally restored. 
It should be noted that in infinite volume limit, the VanderMonde regularization becomes 
essential and, subsequently, Eq.(!67ll is reduced to 



7(11) 
''colorless 



i/3,V) 



yill) 
''colorless 



(/3,^;^o) 



(74) 



where 



y{II) 

Zj 



colorless 



{13, V- a) 



n 



Nc-l 
n=l 



n\ 



exp 



/3C^*(/3,V;a 



/iVe(27r)i(^-i) 



/3 



)(2) 



K^ (/3, V- a) + 2NM,'> (/3, V) 



(2) 



(ivg-i) 

Ireg j 



77'(75) 



where 'jreg = T- and a^ is the lattice space size. Nonetheless, the regularization procedure 
is not required for finite colorless quark and gluon bag. 

The order parameter for GW-like phase transition is the temperature (and the chemical 
potentials). The point of the phase transition, namely. Tow is determined by the continuity 
of the partition function from the low-lying energy solution to the high-lying one and this 
condition is satisfied when both solutions match each others as follows, 



7(1) 
^colorless 



iP,V) 



7(11) 
''colorless 



iP,V) 



(^=I3gw 



^ z 



(I) 

colorless 



(/3Giy,V^;a(^)) 



7(11) 
''colorless 



(/3Giy,V^;cT(^')). 



(76) 



The values of a^^^ and a^^^^ are determined using Eqs. (!53l) and (ITOll . respectively. The 
both solutions (I) and (II) are asymptotic solutions for the low and high temperatures. 
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respectively. This implies that the solution (I)'s partition function is extrapolated to the 
solution (II) when the threshold of GW-like point is reached. Beyond that point, the solution 
(I) deviates significantly from the exact numerical one and turns to be no longer correct. 
Fortunately, the solution (II) provides a clue whereabouts GW-like point threshold and its 
ultimate point and also their interpolation range (i.e. the interval between the threshold and 
ultimate point). The Helmholtz free energy of solution (II) has a hidden valley. The validity 
of the asymptotic solution (II) is maintained whenever the energy climbs the hidden valley 
and reaches the same level of its virtual top that appears at lower temperature. This point 
is ultimate point of the extended interval of GW-like point. Beyond the ultimate point, 
solution (II) matches the exact one precisely. When solution (II) intersects solution (I), the 
threshold of an extended GW-like interval emerges. At GW-threshold, the solution (I) starts 
to deviate significantly and subsequently solution (II) turns to be the correct asymptotic 
solution instead of solution (I). The deviation becomes significant when GW-ultimate point 
is reached. Therefore, it is reasonable to interpolate solution (I) from GW-threshold to the 
asymptotic solution (II) at GW-ultimate point. 

The PNJL-partition function can be solved exactly. The partition function for the color- 
less quark and gluon blob reads 

Zcolorless \P ■, V) = Zcolorless \P i V] Co) ; (77) 

where 

dOi I d92 A^Haar (^i ' ^2) exp [-/3 Vlq^g {[5,V]a,9i,02)\ . (78) 
The invariance Haar measure is given by 

AHaar ("1. 0.) = fg^, n 4 ™' (^) ■ (^j 

Furthermore, chiral mean field, namely, ctq is evaluated by extremizing the partition function 
as follows 

d 1 



loge [Zcolorless (/?, V] Cx)] 



daVP 
Thus Eq.(l80l) is reduced to 

dO^ r d92 /iHaar (^i' ^2) exp [-/3 fi,g, (/3, V- a, ^1, ^2)] 

J —IT 



0. (80) 



X 



da V 
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0. (81) 



Hence by using Eqs. flSOl) and (IHTj) the o"o-chiral mean field is determined by solving the 
following equation, 






da V 



exp[-pn,-^g{/3,V;a,e,,e2)]. (82) 



IV. GW-LIKE POINT AND HAGEDORN STATES 

The asymptotic mass spectral density of states is given by the micro-canonical ensemble. 
The micro-canonical ensemble can be derived from the mixed-grand canonical ensemble of 
a single QG-bag. It is given by the inverse Laplace transform as follows 

Pcolorless (W, \/) ~ -- / d/3 e^ ^ Z.olorless (/?, V^) , (83) 

where W is the energy of QG-bag. In the limit of large W , Eq.( l83|) is evaluated using the 
steepest descent method. The approximation of the steepest descent method fails in the 
limit of small W . This means that it is reasonable to replace the low-lying mass spectral 
density with the discrete mass spectrum of the hadron states while the high-lying mass 
spectral density in the large W hmit is approximated to the bootstrap-like mass spectral 
density for Hagedorn states. Hence, it is more appropriate to replace solution (I) with the 
discrete mass spectrum of hadronic states. Furthermore, it will be shown below that the 
extrapolation of the mass spectral density of solution (I) to Hagedorn states does not lead 
to a deconfinement phase transition to QGP at Hagedorn's temperature. In contrary, the 
mass spectral density for solution (II) leads to a first order phase transition. Therefore, the 
existence of Hagedorn states is interpreted in the term of GW-like phase transition where 
the discrete hadronic mass spectrum turns to the continuous bootstrap-like mass spectrum 
when the hadron's mass exceeds a specific mass threshold (i.e. mu > 2 GeV). In order to 
simplify the calculation drastically, the chiral field is dropped in this section. The density 
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of states for solution (I) is reduced to 

1 r^o+ioo 

1 /"/Jo+Joo 



dPe^'^ziil,^^^W,V) 



271 i 



^o g/3 W^ g-/9VpvdM(/3;'I'o,*o)-/3f^gM(^''^'*0'*o 
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Po^ioo 



y^PVdM['^Q,'^o) rl^O+ioo 



27Ti 



13 W ^^"wsf^o.^oj 



d/3e^'' e 



PQ—ioo 



where W = W + ^^^^^V and 



WpvdM ($0, $o) = 2 loge (l - 6$o$o + 4 (^$1^ + $0) - 3 (<l>o$o) 



(84) 



(85) 



and 



"qqg 



{%,%) = 2Nj 



dxx^ 

27r2 



($0 + 2<l>oe-^) e"^ + e-3^ 
1 + ($0 + $oe-^) e-^ + e-3^ 



($0 + 2$oe-") e-" + e- 



3x 



UJq%% 



(86) 



(87) 



1 + ($0 + ^oe^"^') e-^ + e-3^ 
In the hmit of $0, $0 -^ 0, Eq.([86]) is simphfied to 

lim aggg {<!>0,^q) = "^ f "^ ) • 

<&o,*o^o ^ ^ 81 V 60 / 

Under the assumption of MIT bag model and in the limit of $0) "^o — ^ 0, the extrapolation 
of the mass spectral density (I) is reduced to 



(88) 



where m = W + BV and i?4 ~ 200 — 250 MeV is the bag constant and 



1/2 



33 60 35 

b = l3{i), 
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2V27r \ ^ Iz^ 



(89) 



It is more appropriate to represent the large W limit in the term of asymptotic solution (II). 
Under the assumption of solution (II), Hagedorn's density of states is approximated to 

1 pl3o+ico 



lim pZiLss (W, V) 



2TTi 



d^e^'^ziiH.^ssifi^V). 



(90) 
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In the context of MIT bag model, Eq. (pO|) is reduced to 



where 



Pill) 
b 




iNf + 2-^ NJ-^ 



It is interesting to note that the mass spectral density for solution (I) does not lead to 
a deconfinement phase transition at Hagedorn's temperature while the system with mass 
spectral density (II) undergoes a first order deconfinement phase transition. For a system 
with two flavors (i.e. Nj = 2) and B^'* = 250 MeV, Hagedorn's temperature for the 
deconfinement phase transition is reduced to T^ ~ 608 MeV and 176 MeV for solutions (I) 
and (II), respectively. With smaller bag constant B^^^ = 200 MeV, Hagedorn's temperature 
is reduced to T^ ~ 487 MeV and 141 MeV for solutions (I) and (II), respectively. Hagedorn's 
temperature for the solution (II) is more acceptable than that for solution (I). When the 
exponent a in p(r7i) oc m~"e''''" runs over 5/2 < a < 7/2, Hagedorn matter undergoes 
a higher order phase transition while the system undergoes a first order phase transition 
for 7/2 < a. Therefore, GW-like phase transition is interpreted as an extrapolation of 
the discrete mass spectrum of the conventional hadronic states that are found in the data 
book 311] to Hagedorn states (i.e. super massive hadronic states) that are represented 
by the bootstrap-like models. In this context, the deconfinement phase transition to QGP 
takes place at Hagedorn's temperature. In this sense, GW-like transition is a hadronic 
mechanism that produces meta-stable super-massive hadronic states (known as Hagedorn 
states) below the deconfinement phase transition to QGP. Finally, it should be noted that the 
regularization procedure for VanderMonde's number of states reduces the spectral density 
to 

Piiiliessim) oc m-^e"-. 



oc m 



-(2l^iN^-l) + l)^bm_ ^gg^ 
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Eq.( p3|) demonstrates that jreg may be related to the bag's volume fluctuation. It is reduced 
to 'jreg = 1 for a bag with a sharp surface boundary. The cases 7r.ec, < 1 and •jreg > 1 
correspond to the expanding (dilute) and squeezing (compressed) bags, respectively. The 
case jreg < 1 is related to the bag with an extended surface boundary. The exponent a is 
reduced to | and | for 'j^eg = 1 and 7^63 = \, respectively. 

V. DISCUSSION AND CONCLUSION 

We have considered the canonical ensemble for colorless quark and gluon blob. The col- 
orless quark and gluon blobs emerge as meta-stable Hagedorn states in the relativistic heavy 
ion collisions. These colorless states (i.e. Hagedorn states) significantly enrich the deconfine- 
ment phase transition diagram. Their production signature may mix and be confused with 
QGP. In order to make the discussion simple, at first we neglect the effect of chiral field and 
simply assume massless 2-fiavors in order to simplify the analysis of GW-like phase tran- 
sition. The low-lying energy solution, namely, solution (I), is determined by non-Gaussian 
stationary point method for Polyakov loop parameters ($, $) as defined by Eq. (l33|l . The 
high- lying energy solution, namely, solution (II), is determined by the Gaussian saddle point 
approximation. The assumption is that the non-Gaussian stationary points of solution (I) 
turn to Gaussian saddle points in solution (II). The solution (II) is introduced by Eq.f l51|) . 
Furthermore, the exact numerical solution is considered by evaluating the exact numerical 
integration over the fundamental gauge fields 61 and 62 with the invariance Haar measure 
which is given by Eq. fl271) . 

Fig. [1] depicts the quantity ^ logg Z colorless (T, V) which represents the negative Helmholtz 
free energy density vs T for quark and gluon blob with various volumes R = 0.57, R = 
0.71, R = 0.82 and R = 0.90 fm. It is shown that solution (I) matches exact numerical 
solution below GW-threshold temperature and then it deviates from the exact one when the 
temperature reaches and exceeds GW-threshold point while solution (II) converges to the 
exact numerical solution as temperature approaches GW-ultimate point until it fits precisely 
the exact one as temperature exceeds that point. Therefore, solution (I) is the correct 
asymptotic solution for any temperature below GW-threshold point while solution (II) is the 
correct asymptotic solution for any temperatures above GW-ultimate point. Furthermore, 
it seems that neither solution (I) nor solution (II) fits correctly the exact numerical solution 
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over an extended GW-like point domain which covers the interval between GW-threshold 
and ultimate points. Evidently, the interpolation of both solutions (I) and (II) over the 
interval between GW-threshold and ultimate points fits the exact numerical solution. This 
makes a smooth transition from solution (I) to solution (II) over an extended GW-like point 
interval. The domain between threshold and ultimate points (i.e. over the extended GW-like 
point interval) is reduced to a single point in the limit Nc ^ oo but a finite coupling constant 
g N^. Therefore, the analytical solution is modified smoothly from solution (I) to solution 
(II) over the extended GW-point interval. This implies that low-lying and high-lying mass 
spectra remain in mutual exchange reaction over the extended GW-point interval. The high- 
lying energy solution (II) significantly deviates from the exact solution at temperature below 
GW-threshold and then turns to converge to the exact one as the temperature approaches 
GW-ultimate point and then remains in an excellent match as the temperature increases 
beyond GW-ultimate point. On the other hand, the low-lying energy solution (I) matches the 
exact numerical solution precisely for temperature below GW-threshold and then it starts to 
deviate significantly from the exact one when the temperature exceeds GW-threshold. This 
deviation is significant as temperature increases above GW-ultimate point. The smooth 
modification in the solution's analytic behavior through the extended GW-point interval 
clearly implies that Hagedorn states emerge as meta-stable states over an extended GW- 
point interval with mutual and exchange chemical reaction between the high-lying and low- 
lying hadronic states. The exchange reaction clarifies the difficulty to detect Hagedorn states, 
GW-like transition and the subsequent confusion with the deconfinement phase transition. 
The order parameter $o for solution (I) vs temperature is depicted in Fig. [2] with var- 
ious volumes of colorless quark and gluon bags. The order parameter $o is simply the 
stationary point that projects the color singlet state under the assumption of solution (I). 
The order parameter $o is found very small at low temperatures and this is because of the 
strong confinement. This implies a reduced gluonic component for hadronic states at low 
temperatures since $o smd $o correspond the gluon condensates. Furthermore, when the 
system is heated up, the value of $o increases and approaches its restoration value from 
below but remains $0 < 1- The $o's asymptotic restoration indicates loose confinement 
states or meta-stable bubbles. Furthermore, $0 increases from $0 ~ 0^ to the restoration 
value $0 ~ 1~ within the extended GW-point interval as the asymptotic solution switches 
from (I) to (II). Although, GW-like phase transition from strong coupling to weak cou- 
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pling has been extensively considered in the context of large- A^c limit, it is evident that 
GW-like phase transition persists to exist even in QCD with Nc = 3 but with different 
analytical behavior. The GW-like phase transition in QCD is not a conventional confine- 
ment /deconfinement phase transition but is the Hagedorn's production threshold. This 
can be understood in the term of micro-canonical ensemble and the consideration of gas of 
Hagedorn states. The micro-canonical ensemble of solution (II) is the mass spectral density 
of Hagedorn states where GW-like point corresponds the Hagedorn's mass threshold (i.e. 
^TT-threshoid ~ 2 GcV). Therefore, the present results suggest that GW-like phase transition 
persists to exist in QCD and, subsequently, the Polyakov loop restoration turns to be the 
onset of GW-like phase transition or semi-QGP and the existence of Hagedorn states where 
Hagedorn states are produced in the hadronic phase. This interpretation, definitely, im- 
plies that the Polyakov loop restoration is not the deconfinement 's order parameter as has 
been suggested in some models [l| but the abundant production of (meta-stable) Hagedorn 
states below Hagedorn's temperature. This conclusion is also true for PNJL model where 
the (T-chiral field is considered explicitly and self-consistently in the calculation. Figj3] dis- 
plays the negative Helmholtz free energy vs temperature with various volumes. The general 
situation looks very similar to Fig. [H The solution (I) matches the exact numerical one 
for temperature below GW-threshold. When the temperature exceeds GW-threshold point, 
solution (I) starts to deviate significantly above the exact one. It continues to deviate above 
the exact numerical one as the temperature increases. On the other hand, the high-lying 
energy solution, namely solution (II), has a hidden valley that deviates significantly from the 
exact numerical one at low temperature as far the temperature remains below GW-ultimate 
limit. This valley emerges due to the unphysical oscillatory behavior of the Gaussian sad- 
dle point approximation below GW-like point. Nonetheless, solutions (I) and (II) intersect 
with each others at GW-threshold temperature below GW-ultimate temperature. When 
the temperature increases and reaches GW-ultimate point, solution (II) converges to and 
matches precisely the exact numerical one. Furthermore, as the temperature increases and 
exceeds GW-ultimate point, solution (II) converges to exact numerical one and remains in 
an excellent agreement. Therefore, evidently there is a switch from solution (I) to solution 
(II) within the extended GW-point interval. Nevertheless, the extended GW-point interval 
is ambiguous in heavy ion collisions as far neither solution (I) nor solution (II) fits the exact 
one while their interpolation seems to fit to the exact numerical solution. The importance 

27 



of this inechanism is that it may shed the hght on the existence of (nieta-) Hagedorn states 
and their production as super- massive hadronic states (i.e. ttlh > 2 GeV). The extended 
GW-point interval (i.e. the interval between the threshold and ultimate points) is very sensi- 
tive to the fireball's volume. For instance, the extended GW-point interval is extended from 
Train to T^ax (~ 181 — 291 McV) for bag's radius R = 0.57 fm. The extended GW-point 
interval is significantly reduced and turns to ~ 119 — 184.5 MeV for bag's radius R = 0.90 
fm. Hence, Hagedorn states turn to be of the size of QGP (i.e. R > 0.90 fm) for temperature 
close to the deconfinement one (i.e. T ^ 184.5 MeV). Furthermore, Hagedorn states with 
size R = 0.57 fm are likely to be produced at rather high temperature T^ax = 291 MeV 
while large Hagedorn states are produced at lower temperatures. This unusual behavior 
makes more difficult to detect Hagedorn states as far they emerge as super- massive, glu- 
onic rich and meta-stable states with the size order of QGP. The large (volume and mass) 
Hagedorn states can be confused and mixed with a true deconfinement phase transition's 
signature. The (super-)massive Hagedorn states can be developed as droplets of semi-quark- 
gluon plasma in colorless states. Fig. H] depicts the order parameter $o for solution (I) vs 
temperature with various volumes. The restoration of Polyakov loop $0 likely takes place 
over the extended GW-point interval (i.e. between GW-threshold and ultimate points). 
When T approaches GW-threshold, $0 (and $0) starts its significant restoration process. 
Furthermore, when T exceeds GW-ultimate point $0 turns to be almost restored from below 
(i.e. $0 ^ 1)- This behavior hints that the non-Gaussian stationary point approximation 
fails at temperature above GW-threshold. Subsequently, the non-Gaussian stationary point 
approximation is converted to Gaussian saddle point approximation. 

The effective chiral field G a vsT with various volumes is displayed in Fig. [51 The cr-chiral 
mean field is considered self-consistently in the frame work of solutions (I) and (II) as well 
as the exact numerical solution. The results show clearly that the chiral restoration likely 
takes place within an extend GW-point interval but below GW-ultimate point. Further- 
more, the solution (II) 's chiral restoration takes place before that one for solution (I). This 
implies that chiral restoration likely takes place within an extended GW-point interval above 
GW-threshold point but below GW-ultimate point. The both solutions (I) and (II) fail to 
locate the precise position of chiral restoration. This deficit is understood by realizing that 
extrapolation of solution (I) or (II) is not the correct asymptotic solution over the extended 
GW-point interval. Furthermore, the chiral restoration of the exact numerical solution usu- 
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ally occurs on the right hand side of solutions (I) and (II) but below GW-ultimate point. 
This can be interpreted as exchange reaction and smooth transition between the low-lying 
and high-lying hadronic states over the extended GW-point interval. The chiral restoration 
takes place within an extended GW-point interval but far away from GW-ultimate point 
for small fireball {R ~ 0.57 fm). The restoration point approaches GW-ultimate point from 
below as the fireball size increases. Hence, Hagedorn threshold production takes place within 
an extended GW-point interval. The results demonstrate that the (meta-) Hagedorn states 
are chirally restored above GW-ultimate. 

The results suggest that there is a new class of phase transition in nuclear matter in 
particular in the hadronic sector. The Hadronic phase turns to be smoothly dominated by a 
gas of colorless quark-gluon bags through the extended GW-point interval. This mechanism 
can be understood in the term of Hagedorn states. The size of Hagedorn's bags continue to 
grow up until Hagedorn's temperature is reached. When Hagedorn's temperature is reached, 
the system undergoes a deconfinement phase transition to QGP. The finite volume colorless 
states have been suggested before by Elze, Greiner and Rafelski |l8l-l22l|. This picture has 
been extended to the gas of bags. The results also suggest a possible production of large 
(meta-) colorless quark-gluon droplets of the size order of quark-gluon plasma R > 0.90 
fm at T < 184.5 MeV (in the case of a single droplet analysis). The smaller Hagedorn 
states are produced at much higher temperatures. For instance bags with R ~ 0.57 fm 
are produced at T ~ 291 MeV. This supports that small size hadronic states belong to 
the low-lying hadronic mass spectrum rather than high-lying hadronic mass spectrum. In 
the case of gas of bags, the analysis can be extended using Hagedorn's density of states 
that is derived from the micro-canonical ensemble. This indicates that the nuclear matter 
undergoes smooth transition from low-lying mass spectrum to (meta-) Hagedorn states (or 
even semi-QGP) rather than directly to true deconfined QGP. The deconfinement phase 
transition takes place at Hagedorn's temperature. The production of colorless QG-fireballs 
enriches the nuclear phase transition diagram significantly. This mechanism opens a window 
to produce (meta-) stable colorless super- massive QG-droplets at the size of order of QGP. 
Finally, the signature of deconfined QGP in the heavy ion collisions may be confused and/or 
mixed with the gas of colorless QG-bags or semi-QGP. 
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FIG. 1: (Color online) The negative Helmholtz free energy (density ylog^ Zcoioriess{T,V) for the 
low-lying energy (I) an(i high-lying energy (II) solutions and the exact one (i.e. the phase transition 
solution (I) — )• (II) between GW-threshold and ultimate points) vs temperature. No chiral field is 
considered. The GW-threshold and ultimate points are shown as circles in the T-axis. The results 
for volumes 100, 200, 300 and 400 GeV-^ (i.e. R = 0.57, 0.71, 0.82 and 0.90 fm) are displayed 
respectively in (a), (b), (c) and (d), respectively. 
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FIG. 2: (Color online) The Polyakov loop (triality) parameter <1>, the order parameter of solution 
(I), vs temperature with various bag's volume. The chiral field is not included, (a) V = lOOGeV" . 
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FIG. 3: (Color online) The negative Helmholtz free energy (density y logg Zcoioriess {T, V) for the 
low-lying energy (I) and high-lying energy (II) solutions an(i exact one vs temperature. The a- 
chiral mean field is included simultaneously in the calculation. The GW-threshold and ultimate 
points are shown as circles in the T-axis. (a) V = lOOGeV"^ {R = 0.57fm). (b) V = 200GeV"3 
{R = 0.71fm). (c) V = 300GeV-3 {R = 0.82fm). (d) V = 400GeV~3 {R = 0.90fm). 
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